CRDlKNSWC-HO-1427-01  Iflathematical  Model  for  a  Real  Time  Ship  Maneuvering, 

Stationkeepfng,  and  Seateeplng  TnMng  Simulator 


-  \ 

Carderock  Division 

Naval  Surfeca  Warfare  Center 

Betheada,  Maryland  20084-5000 


AD-A274  514 

I  latWH  tmm  ^ _  •  • 


CRDKNSWC-HD-1427-01  September  1993 
Ship  Hydromechffliics  Department 


Mathematical  Model  for  a  Real  Time  Ship 
Maneuvering,  Stationkeeping,  and 
Seakeeping  Tialning  Simulator 
By 

Paul  J.  Kopp 

OTIC 

lELECTE 
k<iAN04t9S4 


Approved  for  PubUc  Release.  DIstdbution  Unlimited 


xxf^94-00091 

4  tanHi  aim  iiM  aaai  iB»i  mai  Bfi  nai 


1  03  039 


CODE  Oil  DIRECTOR  OF  TECHNOLOGY,  c  J^NS  AND  ASSESSMENT 
12  SHlr*  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

1 .  DTRC  Imports,  ■  formal  aaries,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  re9ardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semiformal  aeries,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3  Technicid  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case-by-case  basis. 


SMUWrrOAMWieATIQNOl'TWtMGHI 


REPOflT  DOCUMENTATION  PAGE 

1  Form  Af^nwad 

1  orye  No.  0704-0188 

la.  ie*0Rr  sEomnY  cuuoHncAncM 

UNCLASSIRED 

tb.  RESnvenvE  MARKmOS 

aa  SECUWTTCUSSIFICATICN  AUTHORITY 

2b.  CBXAS«RCATX)WX»IIWIQIWDI>« 

Approved  w  PutNlc  Release, 

Dtstrfoution  Unlimited 

A  PStraHMMQ  onoANCVmON  REPORT  NUMBSV^ 

CRI:KNSWC-HD-1427-01 

S.  MOMTORMQ  OnOWdZATlON  REPORT  NUMBER 

&  NAME  OF  PBWORMma  ORGMMZAT10N 

Canfercck  Divisian,  Navai  Surf^e 
Warfare  Cer^er 

OFFICE  SVIlSOL 
AfMDictfliM 

1561 

Ta.NAM^OFMONlfonNaoniM«2ATKlN  , 

«c.  AOORESS  (OY.  aaN.  and  ZIP  Coda) 

Tb.  ADDRESS  (Cay,  Stab,  and  2P  Ceda) 

Bethesda,  Mstfyland  20084-5000 

Al  NAME  Of: 

OnOAMZATION 

Coa^ai  Sy^ems  Station 

•MomcEsyMeoL 

ffanUcabK 

3120 

8.  PROCURatlENT  mSTUFiUan’  NUM» 

Ac.  ADDRESS  {Cl^>  SIM*,  and  aP 

ia  SOURCE  OF  FUNDINQ  NUMBBTS 

Panama  City,  Rorida 

PROGRAM 

BSeiTNO. 

63502N 

PADJECT 

NO. 

TASK 

NO. 

WOmiMtTS 

ACCESSONMO. 

DN503466 

a  1  ■  T1TI£  StoMit)  rHwWr  fDan) 


Mathematical  Model  for  a  Real  Time  Ship  Maneuvering,  Stationkeeping,  and  Seakeeping  Training  Simulator 

1Z  P0)SOF^AUTHOn(S) 


Paul  J.  Kopp 


13a.  TYPE  OF  REPORT 

13b.  time  COVERED 

1 A  DATE  OF  Re>onT  (Taar.  Mm«x  Oi«) 

15.  PM3E  COUNT 

Final 

FROM  1/93  TO  7/93 

September  1993 

30 

IS.  SOPPLEMEKTAflV  notation 


17.  OOSATI  COOES 

1 1.  SUBJECT  '^^UIS  ^erttriM  on  fwwM  tf  nooaoMfy  «nd  WarMfy  by  Uodt  nurvibor) 

Maneuvering  Seakeeping 

Stationkeeping  Simulation 

FIELD 

GWDUP 

SUBGROUP 

1 S.  ABSTRACT  (Corbnua  on  lavana  d  naeataaiy  and  WarWy  by  block  lurbw) 

This  report  documents  the  development  of  the  mathematical  model  of  a  combined  maneuvering,  stationkeeping, 
and  seakeeping  simulation  computer  program.  The  model  features  a  mbt  of  physical  based  models  for  the  propulsion 
systems  and  control  surfaces  and  hydrodynamic  coefficient  based  hull  forces.  The  approach  follows  the  concept 
of  the  modular  maneuvering  model.  A  quasi-steady  assumption  is  utilized  for  the  seakeeping  motion  effects.  This 
allows  the  calm  water  maneuvering  rm^ons  to  be  calculated  separately,  and  the  six  degree  of  freedom  lineeu 
response  in  waves  superimposed.  The  effects  of  wind,  current,  and  second  order  mean  drift  forces  are  included. 


20.  nSTRIBUTION/AVMLABIUTY  OF  MSTRACT 

Unlinked 

21  abstract  sbcurry  classification 

UNCLASSIFIED 

22a.  NAME  OF  RESPONSieLE  mOVIOUAL 

Paul  J.  Kopp 

32b.  TREPHONE  (Hvluba  Ana  Codt) 

(301)227-5119 

22c.  office  SYMBOL 

1561 

DO  Femi  1473,  AJN  W 


I 


UNCLASSIFIED 


Contents 

Abjarart  . 1 

Aibi^xtritfhrelnformadoQ . . . 1 

Introthictm  . 1 

Coiid>ixied  Maneuvedi^  Stationke^nng,  and  Scakequng . 3 

Ckxmlinate  Systems . . . 4 

Maneuverii^  Equations  of  Motion  . 4 

External  Forces  . 5 

Curreit  Effects . 6 

HunF(m»Model  . ..7 

Rudder  Force  Modd  . 8 

Pr(q>eller  Force  Model  .  10 

Bow  Thruster  Force  Model  .  12 

Wind  Force  Modd  .  12 

Wave  Drift  Force  Modd  .  13 

Engine  Mscfabery  Modd  .  14 

Seakeeping  Effects .  15 

Summary  .  17 

References .  23 

i’igures 

1.  Conceptual  Modd .  19 

2.  Coordinate  Systans  and  Sign  Conventions .  20 

3.  Rudder  Inflow  Angle  Geometiy  .  21 


Notation 


K 

At 

Cd 

Cd. 


g 

I. 

J 


h 

Jt 

Kq 

Kt 

L 

K 


l^axtt  srea  of  propdlfir  race  avex  rudda 
area  of  rudder 

ciDSisectiooal  area  of  thruster  opeoing 
airfoQ  drag  corfficient 

noQ-dioiensional  wave  drift  momnit  isJMXit  sh^)  z  axis 
rudder  drag  coeffident,  coplanar  to  rudder  ceoteriine 
UHi^iimendoixal  wave  drift  force  in  direction  of  ship  x  axis 
non-dimendotud  wave  drift  force  in  direction  of  ship  y  axis 
airfoil  lift  coeffident 

rudder  lift  coeffident,  normal  to  rudder  anface 
propdler  diameto' 
moment  about  ship  z  axis 
force  in  ship  x  axis  direction 
force  in  diip  y  axis  direction 
gravitational  acceleration 
mass  momoit  of  inertia 
advance  coeffident 
advance  co^dent  based  on  torque 
advance  coeffident  based  on  thrust 
open  water  propeller  torque  coeffidrat 
open  wato’  propeller  thrust  coeffident 
ship  length 

turning  radius  to  ship  center  of  gravity 


iv 


S5(w) 


tun^  nidhis  to  pdiit  OG  the  ship  c<sita^ 
tumisg  n£us  to  tK)i]:it  off  the  ^  ceateriij^ 
tiin^  mSus  to  niddo: 
tiopioass 

ievoliitit»is  per  smmte 
Fevolutbos  jMr  second 
iropdler  |dtd)/diaiBeter  ratio 
efSective  iKusepowo’ (EHP) 
torque 

resistaiice  force 

transfia*  fonction  amplitude 

yaw  rate 

wave  spectra 

thrust 

wave  modal  period 
noth<iiinaiaonal  wave  modal  p«iod 
ship  draft 

surge  vdoctty,  in  ship  coordinate  ^^stem 


DTK:  QUALITY  INaPKinED  3 


surge  vdodty  relative  to  current,  in  ship  coordinate  system 
surge  vdodty  relative  to  wind,  in  ship  coordinate  s;^an 
vdodty,  gnoal 
q)eed  of  advance 

diip  vdodty  magnitude  relative  to  current 
vdodty  of  current 


iootBalon  For 

»TIS  QRAtl  rW 

DIIC  TAB  Q 

fcfiOQOttXlOBti  Q 

Juatlfioatlan _ 

By- _ _ 

Dla^rlbatlao/ 
Svallabiltty  C«nIm 
javail  ~ai^or 

I'  I 


IH'I 


V 


vdodty  at  njd4«^ 

^active  velocity  ov^*  rudda* 


i^lodtyofmnd 


V 
Ve 

V. 

V 

x.y 

x’.y* 

Vit 

y-p 

<i> 

p 

Pc 

PcL 

Pn 

Pnid 

V 
b 


sw^r  velocity,  m  ship  cooniinate  tystem 

sway  vdodtyrdadve  to  curraot,  in  ship  coordinate  tystem 

sway  vdodty  idative  to  wind,  in  ship  coofdinato  syi^em 

squoition  didance  b^eea  propeller  and  fuddff 

rudder  didance  from  ship  carter  of  gravity  (Native  forward) 

thruster  distance  from  ship  cmitor  of  gravity  (Native  torwaid) 

earth  fixed  coordinate  syston  ax^ 

coordinates  of  diip  ceit^  of  grasity  in  earth  fixed  coordinate  systen 
ship  coortfinate  systmn  axes 
translating  coordinate  system  axes 

distance  of  pc^  off  the  ship  centoiine  (positive  to  starboard) 

sq>er8tion  distance  bdw^  propdlms 

blending  function  for  ship  speed 

bloxiing  function  for  relative  heading  to  the  waves 

drift  angle 

drift  angle  rdative  to  current 

drift  angle  at  point  on  the  ship  centerline 

drift  angle  at  point  off  the  ship  cenietiine 

drift  angle  at  rudder 

ramiom  phase  cngle 

tudder  angle 


vi 


£kiw  8xi^e 

infill  fij&dm  phase  angle 
wftvelieii^ 
sigiBScantiwwel^ght 

Aifg^  hesv^  yaw,  roll,  and  i»toh  is  seakeqssg  coondisate 

water 

ccKiqiass  angle  of  direction  of  current  flow 

oon^iass  an^e  of  direction  of  wind 

conqtass  angle  of  direction  of  predononate  waves 

yaw  81^ 

wavefrecpimcy 

encounter  frequi^icy 


I’W^  thrust  wake  fraction 

I-Wq  torque  wake  fraction 

thrust  dechu^on  frctor 


I-t 


Abstract 


this  rept^  doaonents  maihmatiajl  model  of  a  ambined 
mcateuvering,  stationkeeping,  and  xakeepit^  smuiaticm  caa^mter 
^ogtcBiL  The  model  features  a  mix  <f  ph)mad  based  models  for  die 
jmcpul^m  sj^ms  amtrol  surfaces  and  h^dhr^nmnic  co^derU 
based  huU  fmces.  The  approath  folkms  the  cemc^  if  die  nrnduko' 
maneuverir^  model  A  guad-steaefy  assunqftion  is  utilized  for  the 
secdieqnng  nuMon  effects  This  aUaws  die  calm  -water  maneuvering 
modem  to  be  calculated  separably,  and  da  six  degree  for  freedem  Imewr 
response  dt  waves  stgrerittgiosed  The  effects  of  -wind,  current,  endseamd 
iurdm^  mean  drift  farces  are  iraluded 


Administrative  Information 

This  woric  vras  funded  by  the  Coastal  Systons  Station,  Panama  City,  Florida,  und«^  Work 
Request  WX3i004  and  is  identified  at  CDNSWC  by  Job  Order  Number  1-1561-3 14. 

Introduction 

The  Carderodc  Division  of  the  Naval  Surfece  Warfiue  Oatter  (CDNSWC)  was  tasked  by 
the  Coastal  Systnns  Station  (CSS),  Dahlgren  Divirion,  Naval  Surfiice  War&re  Centn-,  in  Panama 
City,  Florida  to  devdop  a  simulator  model  for  ship  mannivering  and  stationkeqnng,  including 
seake^ing  ^ects.  The  rimulator  is  to  be  incorporate  mto  a  team  training  simulator  for  the 
MCM-1  (Avenger)  class  of  mine  hunter  vessels.  Training  scenarios  will  cover  all  aspects  of  the 
MCM  mission;  mine  sweeping,  neutralization,  and  transiting.  As  such,  the  operator  of  the 
simulate  MCM  will  be  require  to  independently  op  '<ri!te  both  the  controllable/rev«'rible  pitch 
propellers,  bow  thiusto',  and  rudders  over  the  full  spee  range  of  the  MCM  in  wind,  waves,  and 
currttit. 


The  baric  physical  requirements  for  the  combine  maneuvering  simulator  were  to  operate 
in  real  time  on  an  Intd  80386/80387  personal  computer  operating  at  25  MHz,  and  be  written  in 
the  FORTRAN  programming  language  using  no  machine  specific  modules.  The  eventual  platform 
for  the  simulator  program  will  be  a  Silicon  Graphics  model  4440  woricstation  opo^ting  at 
40MHz.  Tte  source  code  will  be  translated  by  CSS  into  the  ADA  programming  language  for  final 


»rii!skwi^tihiet<amtnfeCTiystea.lgtl»findffl«a,tteggmlati(mwiBbertffiBii^«tal6Hz 
sa2i^  ntt  20  to  30  miiiisecauis  aSow^  tin»  st^. 

a  ccm^iiter  i^ogram  to  satisfy  slated  reqi^etnoxts,  two  smarting 

poiids  wei«  evahiatod.  Pint,  <me  of  tbe  existk^  matMuverb^  I»x)granis  ftvatialde  could  b«  used 
ttsi  extei^ve  modificatkKts  made  to  it  In  order  to  siQ)p(»t  ti»  zero  speed  ^>edal  case  f(^ 
ttad(»dcee]Ui%  [1^3]*.  Second,  a  new  |^rog^un  could  be  vnittea  fixns  scrsk^  I^her  option 
wcniy  hiwe  to  be  coupled  with  a  wave  induced  ship  motion  calculation  metitodology.  The  second 
option  of  amsw  prc^ram  was  sdected  becau%  intiintic  support  for  zero  q>eed  ai^  revme  and 
tidewiae  mdions  ccHiid  be  bdh  directly  into  the  niod«d.  In  addition,  a  new  program  could  be 
written  in  a  streamlined,  rnmimal  form  for  computational  ^>eed,  vdudt  is  ease:  than  trimming 
down  existing  program  code. 

In  develo{»ng  tins  smulator  model,  there  were  several  otimr  contidoatioi^  aside  from  the 
basic  contractual  ^jedfications.  The  developm«it  of  tlus  Emulator  was  to  be  ti»  first  opportunity 
at  CDNSWC  to  develop  a  maneuvering  modd  based  on  the  modular  manoivering  mode!  com:ept 
The  modular  modd  coiKq>t  diSers  firom  the  traditional  approudi  in  several  ways.  The  traditional 
approadi  is  to  use  hydrodynamic  coeffideots  derive  from  a  Taylor  smes  mqpantion  of  the  forces 
and  moments  acting  on  a  vessd  [4].  The  modular  modeling  reproach  treats  the  componmits  of  the 
ibrc«s  and  momoits  sq)aratdy.  For  example,  the  forces  due  to  the  rudder  are  composed  of  tlw 
basic  forces  from  the  airfoil  in  a  flow,  the  int«'action  of  the  hull  on  the  rudder  (mqjressed  as  flow 
blodcage  and  straightening  effects),  and  the  intnaction  of  the  propeUo’  on  the  rudder  (through 
flow  accderation  and  straightening  due  to  the  propeller  race).  The  forces  on  the  hull  would 
likewise  be  comprised  of  the  bare  bull  forces  and  the  interactions  of  the  rudder  aid  propdler  on 
the  hull.  Implicit  in  the  modular  model  concept  is  the  use  of  physical  based  ttKdds  for  each  of  the 
components  [S], 

The  modular  maneuvorii^  model  offers  several  advantages  over  the  coeffidott  based 
model.  Since  each  component  effecting  maneuvering  performance  is  sqiarate  in  the  modular 
model,  changes  can  be  isolated.  This  includes  changes  to  program  source  code  and  model  as  well 

'KiiAnaoai  b  fcndcrttm  tiiUd  (t  th*  ail  of  repoit 
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as  dbsii^  to  the  dc^gn.  For  ^sample,  the  efi^  of  s  diSbeirt  »ze  or  sl^pe  rudd^  U  isolated 

to  tl»»  rt^dtt’  model  ami  tlm  mteracdoos  o&  oths*  con^onmits.  Hem^e,  the  ImQ  fierce  model  or 
{^(^>dkr  ^cesmdd  remam  mmhaztged.  TMs  is  pardculaiiy  tu^il  is  pr^imioaiy  desgn  \^diere 
differettf  coocqsta  anizmodgated.  in  latm'phas^  of  <^ga,  modb  t^  data  nu^be  a^^dlablefor 
a  ba^iiie  (kagn.  Hie  doagn  can  then  be  modified  with  diffensot  i^opdlers  or  nidders  without 
seMssatily  invalidadz^  the  prevkms  model  test  woric. 

By  usang  modular  manoivering  model  concqit  for  this  simulator,  the  program  could  be 

written  in  a  hi^y  structmred  manner  with  dear  di^dsions  between  force  and  moment  contributors. 
This  fodlitatM  beder  internal  documentation  of  the  smtrce  code  and  allows  eaid^  mainUmasce 
and  future  expansion  of  capabilities. 

This  Tep<xt  documents  the  development  of  the  combined  mai^uvering,  stationkeq>ing,  and 
seakeeping  stnailatioa  mathematical  modds  us^  in  the  ccm^mter  program.  A  separate 
programme  manual  dooinmnts  the  computer  program  source  code  as  implemmited  for  the 
MCM  class*. 

Combined  Maneuvering,  Stationkeeping, 

and  Seakeeping 

In  developing  a  time  domain  model  of  ship  maneuvering  in  waves,  there  are  several 
probiens  to  be  considered.  There  are  frequency  dependent  coefficients  in  the  equations  of  motion, 
unready  ship  speed  and  heading,  effects  of  the  past  history  of  a  body’s  motion  (mwnory  effects), 
and  fird  and  secomt  order  wave  induced  oscillatory  exciting  forces.  Previous  lUtempts  to  modd 
dup  manoivering  in  waves  have  generally  focused  on  the  specific  problons  of  course  keqiing 
ability,  or  ciqisizing  and  broaidung.  Few  complete  modds  that  can  r^resent  g«iaal  man^v^ing 
in  wav^  have  been  devdoped,  with  even  fewer  implementations  of  those  modds  [6],  In  geieral, 
arch  modds  are  Urge,  cumb«a)me,  expensive,  and  not  well  suited  to  prdiminary  design 
applications  ch-  real  time,  man-in-the-loop  dmuiations. 


*  ts  t  UMii&atMa. 


The  if^oadi  takoi  tee  is  to  decouple  the  rdativt^  slow  aanimvaring  motions  &om  the 
rdativety  &st«ar  seai»qring  motions.  Ttus  is  accomplished  by  calcukdi^  the  calm  water 
maoeuvah^  motions  ai^  adding  the  linear  seakeeping  responses.  The  advantage  of  this  method  is 
that  the  seakeeping  motions  can  be  pre>computed  outside  of  the  maimuv&ing  simulation.  Figure  1 
shows  the  concQ>tual  model  that  was  developed  as  a  framework  for  the  devdopmsit  of  the 
compontes  of  the  modd. 

It  is  important  to  note  that  estimation  or  experimeatal  determination  of  the  hydrodynamic 
duractoistics  of  each  conqsonent  of  the  ship  (hull,  mddm’s,  props,  etc.),  it  is  assumed  that  the 
body  is  mo^nng  with  su£Scient  speed  and  has  sufficient  surface  roughness  to  irrsure  that  the  flow  is 
fully  turbulmit.  However,  the  simulation  model  has  the  c^ability  to  operate  at  very  low  speeds 
whm  the  flow  could  be  lamiruu.  No  attempt  is  made  in  this  model  to  account  for  hydrodynamic 
differences  between  laminar  and  turbulent  conditions. 

Coordinate  Systems 

Within  this  simulator  model,  a  multitude  of  coordinate  systems  and  firames  of  reference  are 
employed.  The  standard  coordinate  systems  for  ship  motions  and  maneuvering  are  different  from 
each  other  by  a  180  degree  rotation  about  the  ship's  x  a^ds.  The  bottom  of  Figure  2  shows  the 
orientadon  of  both  coordinate  systems  [4].  In  this  work,  the  standard  maneuvering  coordinate 
system  will  be  used  unless  explicitly  stated.  Figure  2  shows  the  translating  ship  coordinate  system 
(x’.y'),  the  earth  fixed  coordinate  system  (X',Y'),  and  the  ship  reference  coordinate  system  (x,y). 
The  positive  sense  of  the  rudder  and  drift  angles  and  propeller  and  bow  thruster  thrust  are  also 
shown.  The  origin  of  the  ship  coordinate  systems  are  at  the  ship's  center  of  gravity.  Current  and 
wbd  speed  and  direction  and  predominate  wave  direction  are  specified  in  earth  fixed  coordmate 
system  using  compass  heading  angles. 

Maneuvering  Equations  of  Motion 

The  general,  nonlinear  equations  of  motion  for  a  maneuvering  ship  in  the  horizontal  plane 
(without  roll)  are  given  in  Reference  4 , 
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m(u- v\J>)«2Fx 

m(v  +  u4>) «  2Fy 
» IMi 

Hie  aos8  product  terns  on  the  left  dde  of  the  x  and  y  equations  are  tlw  result  of  these  equations 
being  written  in  the  diip  coordinate  system  Rewriting  these  as  a  set  of  first  order  differoitial 
equations, 

v=^-ur  +  2Fy) 

(2) 

where  hasbeenrqilacedbyr. 

Relating  the  ship  coordinate  system  back  to  the  earth  fixed  coordinate  system  yields  another 
set  of  first  order  difiTerential  equations, 

Xpo,  -u  smijj+v  cosiJ> 

Ypc,  “U  cos  ii>-v  simi)  (3). 

Equations  2  and  3,  when  taken  together,  form  a  set  of  five,  first  ordo*  nonlinear  state 
equations.  The  general  form  is  given  by 

X»A{X)+B(U,X)  (4). 

Solution  of  the  time  evolution  of  these  equations  is  performed  by  first  order  Euler  time  step, 

X»*‘  -X“  +  [A(X")  +  B(U",X“)]At  (5). 

The  first  order  Euler  method  is  generally  stable  (for  dynamically  stable  ships)  and  will  not 
propagate  errors  if  the  time  step  size  is  small  compared  to  the  speed  of  ship  manmivering 
motions. 


External  Forces 

If  we  disregard  the  seakeeping  aspects  of  the  problem  for  now,  the  total  external  forces  and 
moments  acting  on  a  maneuvering  ship  are  represented  by  XFs,  XFy,  and  XFg .  These  summations 
can  be  expanded  into  components  for  the  hull,  rudders,  propellers,  thrusters,  wind,  and  waves  [4], 
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(6). 


ZFy  -  Fy^  +  Fy^  +  Fy^^  +  Fy.^  +  Fy,j^  +  Fy,^ 

ZFb  •“  Fqi^  +  ^®»»in 

Ciirrent  Effects 

HuU  fiovM  are  due  to  the  rotation  and  translation  of  the  v^sel  throu^  the  watn*.  In  order 
to  inchtde  the  effects  of  a  curmit,  the  rdative  velocity  betweor  tlw  ship  and  the  current  must  be 
calculated. 


Uc“VC0Sp  +  V«Bw»  C0S(*«itrt-11>) 

Ve  =  -VsinP  -  Vcwwot  an  (♦«««*  -'!>)  (7) 

where  the  ship  vdocity  and  drift  angle  are 

p  -  ijj-arctan^  (8) 

and  the  drift  angle  is  considered  zero  if  both  the  Xp„  and  vdocities  are  zero.  The  vdodty  and 

drift  of  the  ship  rdative  to  the  current  is  therefore 

Vc-Vu?+v| 

Pc-arctan~  (9). 

Using  the  relative  vdocity  through  the  water  requires  that  the  surge  and  sway  equations  in 
equation  2  be  rewritten  as 

Ue  =  s(Ver  +  2F*) 

Vc  =  s(-Ucr  +  2Fy)  (10). 

Equation  3  must  then  be  corrected  for  the  current  and  written  as 

Xpoi  -usinijj+vcostl> 

Ypoi  -ucos  q»-vsinq>  (11) 

^ere 

U-Ue-Veuw»itCOS(\j)  -  4wiwat) 

V"V6+VeuffiBaiSin('IJ)  — (1^)- 
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Him  Force  Model 


As  stated  in  the  introduction,  the  traditional  approach  to  modeling  hull  forces  has  been  to 
the  lineaiity  of  foic«t  and  expand  them  with  a  Tajdor  series  eiqpandon.  This  is  still  the 
most  straight  forward  ^proadi  to  developing  hull  forces  as  a  function  of  surge,  sway,  and  yaw 
motions.  Here  we  ^  iu>t  indude  the  effects  of  rudder,  props,  or  bow  thrusters  in  the  expansion 
by  conridering  the  bare  hull  only. 


If  the  hull  forces  are  considered  to  be  functions  of  nirge,  sway,  and  yaw  motion  only ,  we 
may  write 


'  air”  « 

T?  A  1  i  ^y  •  i 

t  ■3^1+ 


(13) 


^ere  the  higher  order  tcnm  have  been  omitted.  The  derivatives,  also  known  as  hydrodynamic  or 
maneuvering  co^dents,  are  usually  written  in  a  shorthand  notation  such  that 
~  -Xv,  ^  -Yv,  and  ^  -Nv,  for  example.  The  acceleration  terms  Yv,  Y  r,  N  v,  and  Nj 
basically  represent  low  frequency  added  mass  effects  [4]. 


Hydrodynamic  coeffid«:ts  can  be  obtained  firom  crq>tive  modd  tests,  estimation,  or  system 
identification  fi’om  full  scale  trials  or  fi'ee  running  modd  experiments.  Each  of  these  methods  has 
their  own  merits  and  drawbacks.  A  discussion  of  these  is  beyond  the  scope  of  this  document  and 
the  reado’  is  referred  to  the  references  [4, 7,8,9]. 


The  normal  practice  is  to  express  the  hydrodynamic  coefficients  in  non-dimensional  form, 
for  exanq)ie 


(14). 


Here,  the  vdodty  in  the  denominator  is  the  relative  velocity  through  the  water.  When  a 
non-dimensional  hydrodynamic  coefficient  is  used  to  calculate  the  non-dimensional  force  or 
moment,  it  implies  that  the  assodated  variable  (v  in  the  example  case  of  equation  14)  to  be 
multiplied  is  also  in  its  non-dimensional  form.  Hence,  the  dimensional  force  is  obtained  from 


7 


gumming  the  products  of  the  iK)ii>^imenaongl  co^deots  and  associated  son-dimmisional 
'rariables,  and  then  muh^tying  by  the  factor  (L^  for  nmments) 


\^^thm  the  series  esqpanrion  for  the  x  force,  a  term  has  beoi  added  to  account  for  the  fauQ 
redstance.  The  effective  power  (EHP)  curve  is  rqjresntted  as  a  polynomial  curve  from  ^ch  the 
hull  resistance  is  obtained. 


Pe(V)*  V3  [Ci V  +  C2V*  +  C3  ] 


Rt- 


s5ori 


(15). 


The  polynomial  coeffidohs  are  obtained  from  curve  fits  of  upright  straight  line  (zero  drift  angle) 
redstance  tests  [10].  The  ship  speed  used  in  tlus  inqilementation  will  be  x  conqK}nntt  of  the 
rdative  speed  through  the  wato^  (uj.  In  the  special  case  of  zero  ship  speed,  the  redstance  is  zero. 


Rudder  Force  Model 


The  rudder  is  the  primary  control  surface  by  which  directional  control  is  maintained.  In 
order  to  det^mine  the  forces  acting  on  the  rudder,  the  actual  flow  velocity  and  direction  at  the 
rudder  must  be  determined  from  the  general  motion  of  the  ship.  Figure  3  shows  a  schematic 
representation  of  a  ship  in  a  turn.  If  the  velocity  components  of  the  ship  center  of  gravity  (eg)  and 
the  yaw  rale  are  known  then  the  drift  angle,  Uiming  radius,  and  inflow  angle  at  the  rudder  is 
calculated  from 


Lc^=Ve/r _ 

Ln«d  -sign(r)./(XroicosPc)*  +  (L,,-x«dsin  pc)’ 

-  p-dgn(r)arctan(j;^^) 

V rod  “  IrLnid  I 

a-Pwd  +  6-e  (16). 

In  the  case  of  r*0  and  V»*  0,  then  the  turning  radius  is  infinite  and  pmd  ••  Pe  and  The 

rudder  angle  is  givm  by  the  commanded  rudder  angle  up  to  the  point  that  a  set  maximum  rudder 

rate  is  achieved. 


When  the  point  is  not  on  the  ship  centerline,  an  additional  correction  is  required.  If  the 
turning  radius  and  drift  angle  of  a  point  on  the  centeriine  is  L(^  and  Pcl  >  thon  tho  turning  radius 
and  drift  angle  of  a  point  off  the  centerline  is 
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Lpt  -sign{LcL)^(yptSmpcL)  +  (Lcl  -  YpiCOs  pea.) 

Ppt  -  PcL+sign(LcL)arctan  (17). 

The  correctioin,  e  in  equation  16,  is  a  flow  straightening  coeffldent  inhich  rq>resaits  the 
blodkage  of  flw  hull  to  the  flow.  It  has  been  seoi  that  the  flow  straightening  is  primarily  dq^endent 
on  the  geoQtftric  drift  angle  at  the  rudd^,  Pnid  [1 1].  A  rdation  derived  ftom  the  data  pr^ented  in 
Refo'^ace  11  is  enqiloyed, 

t  -Co+Cipnid+C2pi^+C3p^  (18) 

where  it  is  assumed  that  e  •  0  if  jPntji  ^  90  degr^. 

Tha«  is  an  additional  intoaction  ^ect  on  the  effective  flow  vdodty  when  the  propdlo* 
race  is  directed  ova*  the  ruddo-  [12], 

VL  -  ^(1  -  WT)V„«,cosprud  +kv.]’  +  ^1  -  WT)*ViLiCOS*p„^  (19) 

where 

V.  -  -( 1  -  WT)V«dCOS  P«d  +  J(1  -  Wt)^V„mjCOS^Pto|  +  IKto^D^ 
k-0.5+0.26527tanh  ( 1.27751)  +  0. 17533  tanh(2.555i) 

|-2Xp<^/D  (20). 

Some  portion  of  the  propeBer  race  will  be  directed  ovw  the  rudder  when;  u,>0  and  p/d>0 ,  and 
possibly  u^<0  and  p/d>0. 

The  forces  acting  on  the  rudder  are  given  by  the  lift  and  drag  coefiSdents.  As  the  rudder  will 
operate  in  a  diversity  of  flow  conditions,  lift  and  drag  characteristics  are  needed  for  forward, 
reverse,  and  sideways  flow  [13].  Like  the  hull  hydrodynamic  coefficients,  rudder  lift  and  drag 
coeffidents  are  generally  detomined  in  a  fully  turbulent  flow  condition.  These  co^dents  will  be 
used  at  all  tunes  in  the  emulator  model  even  if  the  flow  speed  and  geometry  at  dther  port  or 
starboard  ruddo'  would  in  fact  be  indicative  of  laminar  flow. 

Lift  and  drag  are  defined  respectively  as  parallel  and  normal  to  the  flow.  In  order  to  obtain 
the  forces  acting  on  the  huU,  a  transformation  is  performed  fi'om  the  flow  reference  frame  to  the 
rudder  refermce  frame  and  finally  to  the  moving  ship's  coordinate  system.  This  is  done  using 
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Cl— '^LCOSa+CD^tt 

-Cpcosa-CL^a  (21) 

and 

F»  «  (Cl^cosS  -  CD««.sm6)(ipvLAf ) 

Fy  -  (-Ci>^co85-Cl»h,^6)(jP^A,) 

Fo  *Fy*Xnid  C^)' 

Thou^  die  ludders  OS  the  MCM  operate  together,  each  expoi«u«  a  difiSsr^  inflow 

^n^ry  aiui  potentially  (hfiferoit  propeUcs^  races.  Therefore,  the  forc^  and  moments  should  be 
calculated  &»-  eai^  rudder  s^aiatdy  and  aunmed  together.  The  machineiy  modd  for  the  steering 
gear  in^ses  liniits  <mi  the  extr^n^  of  position  and  the  maximum  rate  of  change  of  portion. 

Propeller  Force  Model 

The  propdlers  are  the  primary  means  for  the  ship  operator  to  control  the  ^>eed  of  his 
vessel.  They  also  pro^e  some  directional  control  whra  the  two  propdUers  are  op«ated 
individually  at  different  speeds  and/or  pitdi  settings. 

Since  the  MCM  has  controllable  reversible  pitch  (CB?)  propellers,  the  foil  four  quadrant 
open  water  propdler  performance  at  all  pitch  settings  must  be  modeled.  The  open  water  propdler 
data  is  rqiresented  as  curves  of  thrust  and  torque  coefficients  ^  and  Kq  respectivdy)  over  a 
range  of  n^ative  and  positive  advance  coefficient,  J  [14],  These  are  defined  as 

(23) 

where  V,  is  the  ^leed  of  advance  through  the  water,  n  is  the  propeller  speed,  D  is  the  propeUn- 
diameter,  T  is  the  propeller  thrust,  and  Q  is  the  propeller  torque.  For  the  purposes  of  this  model, 
the  ^)eed  of  advance  is 

V.-V{1-wt)  (24) 

where  (l^Wy)  is  the  thrust  wake  fraction  and  V  is  the  x  compon^it  of  the  velocity  of  a  point  at  the 
propeller  hub  (refer  to  the  discussion  in  the  pre'rious  rudder  force  section).  ITie  thrust  wake 
fraction  rqiresoits  a  blockage  effect  of  the  hull  on  the  fiow  to  the  propeller  and  is  usually 
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detennined  by  powenag  esqperunoits.  In  straight  line  motion  (r^o  drift  an^e),  the  wake  fta^on 
is  the  n^asured  value  from  expoimoits  [10],  Howev^,  as  the  ship  maneuvers  and  operate  at 
non-zeio  drift  ang^e^  there  will  be  less  blockage  by  the  hull  on  the  flow  to  tlm  propello^.  For 
inflow  angles  kss  than  90  d^rees,  a  cosine  ^uared  correction  is  us«i, 

(1 -wt)  “  an  PjBop+(l  ”  WT)(pj„^)Cos*Pjrep  (25) 

and  (l-Wj)^!  for  angles  greater  thmi  90  d^ees.  Here  ^ftop  is  the  drift  ai^e  of  a  point  at  the 
propdltf  imb.  Siime  the  propeQos  are  opeeted  independently  from  each  other,  tlu»  V  in  equation 
24  and  in  equation  25  should  be  calculated  ft>r  both  port  and  starboard  prop^en. 

A  polyncmnal  curve  fit  for  each  Kj.  and  Kq  curve  is  made  for  each  pttch/diamrt<»^  (P^) 
value  and  both  positive  and  negative  J.  Reference  [IS]  g^es  twdfth  ordsr  polynondal  coeSSdents 
for  the  MCM  dedgn  propeller.  In  order  to  obtain  the  operating  values  for  either  port  or  starboard 
propdler,  the  polynomials  are  evaluated  at  the  ^propriate  value  of  advance  co^dent.  Linear 
int^polation  is  used  between  curves  to  the  actual  propeller  pitdi/diameter  setting.  The  migine 
modd  (to  be  discussed  later)  determines  the  rpm  to  be  used  in  calculating  J  for  the  port  and 
starboard  props. 

Having  the  total  thrust  from  the  K^.  values,  it  must  be  decreased  by  using  the  thrust 
deduction  fimtor.  Thrust  deduction  is  corrected  for  non-zo'o  inflow  acgl^  in  a  siinilar  manner  to 
that  used  on  the  wake  Section  in  equation  25.  Propeller  thrust  acts  only  along  the  longitudinal 
dire(^on  of  the  ship  hufl.  No  oblique  flow  efiect  are  included  except  for  the  wake  fiaction  and 
thrust  deduction  fitotor  corrections.  As  a  result,  and  by  assumption,  no  side  forces  are  exoted. 

The  forces  and  moments  delivered  by  the  propellers  to  the  ship  are  given  by, 

Fx  -KT(l-t)(ipn2D") 

Fy-0.0 

Fx-±5:fFx  (26). 

For  the  yaw  moment,  the  positive  sign  is  taken  for  the  port  propellm’  and  the  n^ative  for  the 
starboard  propdler. 
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Bow  Thruster  Force  Model 


XliA  bow  tbnuter  is  an  auxiliary  device  for  providing  directional  control.  Tl»  MCM  class  of 
vessdb  use  tin  Omni  Thru^™  bow  thrusts.  This  device  is  difSsent  fiom  the  more  convaitiona} 
tunnel  type  thrusters.  The  Omni  Thruster^  uses  a  sea  ch^  in  the  keel  for  a  water  inlet  aixl  then 
^lits  the  £bw  to  eitho*  a  port  or  starboard  side  outlet.  For  the  purposes  of  this  smulator,  a 
relative  simple  modd  wiQ  be  used. 


The  operator  will  issue  a  thruster  command  of  bdween  -^lOO^  and  -100%  thrust  A 
second  ordo^polyiminial  is  us«i  to  transform  the  command  setting  to  the  full  scale  thrust  Tt» 
data  needed  to  perform  the  curve  fit  is  usually  provided  by  the  thrusto-  mamifiacture  and  should 
assume  no  forward  speed.  Forward  speed  degradation  effects  have  bear  studied  by  Chislett  and 
Bjorheden  [16]  and  McCrei^  has  performed  curve  fits  of  the  published  data  [12].  The  thruste^ 
forces  can  be  d^ermined  fiom 


F*  -  0.0 

Fy  -■Itexp(-13.3m^)  +  0.627396m-0.385772m»  +0.124873m3] 

F„  -  [T-KI  -  f)iVy-r)]x^  (27) 


where 


_ 


(28) 


and  T  is  the  thruster  thrust  with  no  ship  velocity.  It  is  assumed  that  m=0  if  T=0.  For  a  ship 
moving  ait,  it  is  assunsed  that  th^e  is  negligible  degradation.  This  is  represented  here  by  m=0  for 
iPei  >  90  d^rees.  Hie  vdodty  and  drift  angle  at  the  thruster  should  also  be  calculated  for  the 
ship  in  a  turn  using  the  method  outlined  in  equation  16. 


Wind  Force  Model 

In  a  manner  amilar  to  the  current  and  relative  speed  through  the  water,  a  transform  between 
the  wind  speed  and  direction  in  the  earth  fixed  reference  system  and  the  ship  refisroice  system  is 
givoiby 
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(29). 


Uw»VcOS§  +  C0S(4wbii-^} 

Vw  » -Vsinp  -  -•>!>) 

Vw  -  +vi 
Pw^arctan^ 

The  fwi^  &ice  tmi  monuot  coeSBdeoU  are  d^ned  as  tt»  sde  force  and  momeot 

(nos-dmwiiK<»adized)actisgontheship  wboi  the  wind  b  blowing  across  tibe  pent  twain.  The 
sir;^  force  and  nuMneoico^deats  are  defined  wind  is  trowing  across  ti»  bow.  The 

total  force  ami  mennent  cmfofinition  fiom  the  wind  is  a  blended  value  of  the  two  sets  of  wind 
force  ami  im>mait  codEdeots  [12], 

F*  -  |pMrVi(Area,C*^cospw) 

Fy  -  |p«Vi(AresyCy^sm^w) 

Fb- jPiirV^(Area*Cax^cos2Pw+AreayCyx^an2pw)  (30). 

The  negative  agns  are  due  to  du  wind  drag  opposing  forward  modoa 

Wave  Drift  Force  Model 

The  wsvt  drift  forces  being  addressed  here  are  the  mean  value,  second  order  wave  drift 
forces.  This  is  an  acctq^le  nuMld  for  use  in  real  time  man'in-the-loop  siisulation  where 
calculation  of  the  fiill  time  varying  second  ordm  forces  would  be  burdoisome.  The  mean  wave 
drift  forces  are  rqnvsented  by 

F*  -  pgB^l/3CD,(T’o,<^  wwvs) 

Fy  “  PSi^l/jCoyCT  «,  fwaw) 

FB-pgL2^IflCD,{T’o.9v«v.)  (31) 

where 

T’o-To/irr:  (32) 

is  ti»  noDHiimaidonal  wave  modal  period  [12]. 

McCrei^  aiui  Jones  have  investigated  the  wave  drift  forces  on  an  MCM  modd  in 
uni-direcdonal  r^;ular  seas  [17].  Since  then,  McCreight  has  developed  an  ui^ublished 
intepolation  method  for  the  original  model  test  data.  The  non-dimmisional  wave  drift  force 
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codicki^  used  m  ll»  siiBiIatkm  cos^nite*  prognun  are  d^ssiiied  by 
data. 


Engine  Machinery  Model 

The  USS  AVEKC^  O^Otd-l)  is  powered  by  four  Wauke^  LN  1616  DSIN  diesd 
eo^mea  geared  to  two  TrassanrafkaoDdaval  reductfon  ^ars.  There  are  also  tnnn  !i^t  load 

mottos  foat  can  be  UMd.  This  amulator  &K>dd[  will  assume  that  tvro  die^  (bivix^ 

^laft  ire  r^sesented  by  cnie  ^^e  large  en^ne  driving  its  own  ^aft  with  no  ooss  ccHjpling 
between  ihafis.  No  modd  for  the  Kght  load  motors  is  included. 

Abaskmoddfin-adiesd  engine  is  a  constant  torque  model  [18].  Within  the 
implementatiott  offoectmatant  torque  modd,  oiguus  throttle  settings  are  givQi  as  0  to  100%  such 
that  100%  throttle  produces  naaxtmumen^e  torque.  limits  are  placed  on  the  rate  of  change  in 
produced  en^ne  torque.  In  addition,  the  shaft  speed  is  limited  to  anuudmum  value  as  is  the 
rotatiemd  acceieation  of  the  diafi.  The  rate  of  dumge  of  propeller  pitch  is  also  limited  to  a  set 
manmiun  value  [19]. 

Fm*  a  omstant  tenque,  the  shaft  speed  is  not  specified  and  nmst  be  badcsolved  unng  kimwn 
diaracteristics  of  the  inopdler  and  the  flow  conditions  at  tte  propdlo’.  Since  this  nK>del  must 
poform  statioQkeejnng  as  wdl  as  geno-al  maneuvering,  than  are  3ev«d  spedal  cases  that  nmst 
be  accounted  for.  In  addition,  since  reversible  propellers  are  used,  the  shaft  speed  should  be 
prevented  fiom  becoming  n^ative.  The  trivial  case  of  zm'o  toique  (corresponding  to  s^o  throttle 
setting)  {Hoduoes  zero  shaft  speed. 

The  bollaid  puQ  cemdition  exists  when  tbm'e  is  a  iK>n>zaro  torque  and  zicro  ship  speed.  H&e, 
the  advai»e  coefficient  J  is  also  zero.  The  values  of  torque  and  thrust  coeffidents  (Kq  aaiid  K^)  are 
finind  directly  fiom  the  opoi  water  propeller  characteristics  at  a  given  propeller  pitdh  setting, 
interpolating  between  discrete  pitch  values  as  necessary.  The  shaft  speed  is  then  givei  by, 

n-j(y(pK«D’)  (33). 
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final  csdsUi  wi^both  the  torque  and  sMp  ^eed  are  iK»i>zero.  Inthis  ca&e  the 

&ctor, 

^^{pD’Vi)  (34) 

Is  cooqHited.  Hus  is  isd(^)eQden£  of  shaft  ^>eed.  Knov^  the  opes  wate'  propdler  diaiacteristics 
as  a  pdiynosiyal  fiaictios  advaiKe  coefiBdeots,  we  can  write  with  use  equarion  34, 

0  ->co  +CtJQ +C2Jq^1  -Q/^pD^V*^  ^  +C3JQ+ ...  +ctjJ^  (35). 


Ti»  real  root  (or  roots)  is  torque  advance  co^cient  at  the  current  condition  of  shi^ 
shaft  t<Mn^,  and  pit^Mfto'lutch.  The  thrust  advance  co^doit  is  given  from. 


Jt“Jq 


t-wr 

1-W15 


(36). 


The  |»qpdler  opes  water  thrust  coefiSdent  can  be  obt&iiuKi  from  the  known  J.,.  ^nd  the  polynomial 
curve  frtrKf.  The  draft  ipeed  is  given  1^, 


(37) 


<»-  equation  33  if  If  the  shaft  speed  or  rate  of  change  of  shaft  speed  is  ftmnd  to  nci  reed  pre^ 

limits,  ttoi  it  is  apjffopriatdy  down  graded  to  M  within  the  limits.  Since  J7  and  JqW^c  ftmnd 
from  relations  indqKsident  of  shaft  spe^  those  values  will  not  change,  hmice  Kg  wiT  'lot  diange 
either.  However,  the  ei^iiM  torque  will  duinge,  and  its  new  value  is  computed  from, 

Q-KQ(pn2D’)  (38). 


Seakeeping  Effects 


The  sx  d^ree  of  freedom  responses  are  governed  by  frequor<c>'  ^icTnain  transfo'  fiinctions, 
miraiiaitaH  iiflmg  a  liitsw  Ship  motions  computo'  program  [20].  Ti'}  ^  histories  for  tl»  six  degrees 
of  freedmn  are  gesMrated  from  the  amplitude  and  phase  information  contained  witinn  the  transfer 
fimetioos, 

nKt)  -  2  [{Ra*  •  Ck)cos(ii>^t  +  Yk  ^  Pf»)]  (39). 

k 


IS 


11^  gives  time  t  iidsae  the  summstioii  is  over  di^nete  £«qiietK^  and 

£r  is  the  traa£^fii&ctioii  amplitude  and  phase  at  the  w«ve  fieqiuKKiy,  <a«  is  the  frequency  of 
etic(»£mcr  attheii^vveve&^ueBcy,  andtk  is  the  wave  hdght  at  the  wave  frequency  de&u^ 

IQT 

A  uni&nn^  dh^rSnited  random  phase  angle,  Yk,  is  indudni.  Equation  39  8ppl}<»  specifrcalfy  to 

iong  orated  mss.  Short  oested  sea  re^ot^s  require  aii  additional  simmation  ovo-  mve 
direetkm  [21], 

There  are  576  individual  six  d^;ree  of  freedom  time  histemes  used.  Th^  are  cmnpd^  of 
ax  dnp  ^>eed8  (0  to  1 0  knots  in  2  knot  h^mnmhsX  twenty  fmir  leadings  (0  to  345  d^ees  in  1 S 
d^ree  imnements),  and  frnir  sigmScant  wave  hoght/modai  peiod  condnnafrcHis;  1  frx^  and  7 
seconds  (Sea  State  2X  3  feet  and  7  seconds  (Sea  State  3),  6.2  feet  and  9  seconds  (Sea  State  4), 
and  10.7  feet  and  9  secoiuis  (Sea  State  5).  For  theMCM,  each  tinK  lustoiy  is  10  mimites  in  length 
at  4  saiQ{des  per  seconds  and  short  crested  seas  are  asauned.  Each  sm  of  time  histones  are 
groiqied  in  data  fries  acemdiog  to  sea  state. 

Access  to  die  tune  history  data,  by  the  maneuvering  simulator,  is  p«formed  usng  a  data 
table  lot^Hip  vdth  liimar  intepolation  betwe^  beading  and  speed  as  needed.  However,  h  stould 
be  understood  that  linear  ship  motions  theory  assumes  constant  mean  heading  and  constant  mean 
ship  ^)eed.  Both  the  transfer  functions  and  generated  time  histories  are  sensitive  to  ^>eed  and 
heading  ^uiarimis.  Hence,  an  interpoiated  time  lustory  for  a  given  mode  of  motion  may  iX)t 
itfcessaiily  be  the  state  as  a  time  lustory  doived  from  a  transfo-  fui^on  at  that  identical  speed 
and  heading. 

Interpolation  between  q>eed  and  heading  is  performed  using  linear  blending  fru^mis. 

Tirae  fiinctkms  are  dqjeruieot  only  on  speed  or  beading  and  can  ther^ore  be  used  for  all  modes 

of  motion  withmit  recom|xitmg.  For  example,  the  ship  speed  ami  relative  heading  to  the  waves  are 
given  at  disoete  values  and  Xi  ,...,Xn  respcctivdy.  The  blending  funedens  fr;r  ship  speed 
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(0)  ttodl  to  speed  V  aod  X  betwecsi  the  aid 

j^)eed  «a4i*' tod  liNi^i^  tre  ths^ 

<I>(i)-l-.®(i+i) 

'PC) - 1-- ^0+1)  (41). 

Meipolstiag  <t  time  t  takes  place  betweoi  time  hi^ory  values  at  ^leedsi  and  H-l  and  heading] 
mdj't'l, 

rit,  V.  X)  -  [n(t,  Vj^i.  Vi)  •  ’^0+1)  +  n(t.  Vj,j.  X^)  •  ^0)3  •  4>(H1)  + 

Vi.  Vi  )  •  'P(i+1) + n(t.  Vi,  •  ¥0)3  •  ^0  (42). 

Summary 

A  desciiptimt  of  the  mathematical  model  used  in  develops^  a  asmlator  for 
maneuverii^  rtatkadceqwaft  and  seakequng  has  been  presented.  The  modd  is  based  <m  the 
c«ic^  of  a  modular  sb^  maneuvoing  modd.  Hus  treats  tlm  ship  huH^  mddeis,  propdiers,  and 
bow  thruster  as  mdivkhial  eompotusiits.  modded  separatd/  usii^  metlmds  ap^cpriate  fin*  eadi 
ctmipmieDt  with  imeraoitm  effects  between  componosts  accounted  for. 

Future  additjcna  and  modification  to  the  modd  are  anticipated.  The  list  indudes  ti» 
additim)  of  red!  in  the  equadons  of  motion,  the  use  voticai  axis  propdiers,  towing  forces,  shallow 
water  eSsets,  fo^diip  interaction  effects,  ship-foore  impaction  effects,  time  varying  second 
<^der  wave  drift  finces,  and  multiple  imgine  modds.  Anotho*  important  nK}dification  wmild  be  a 
more  riskmd  rqrfnoach  to  the  siq>erpoation  of  reakeq)ing  motions  on  the  cahn  water 
maneuvmiog  motimis.  Hie  conqnita'  program  \diich  in^loiwnts  the  shnulator  nmdd  might  be 
gh«n  sitomatic  cootrd  systmns  for  an  auto  pilot  as  wdl  as  dynamic  poshionkeqiing.  New  front 
Olds  to  the  Hnnilator  will  also  be  devdop^  allowing  for  a  fully  graphical,  roan-in>tl^loop 
sumiiatuHL 
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Figure  1  -  Conceptual  Modd 
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Figure  3  -  Rudder  Inflow  Angle  Geometiy 
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